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Abstract - We discuss the effects of particle exchange through open boundaries and the in¬ 
duced drive on the phase structure and condensation phenomena of a stochastic transport process 
with tunable short-range interactions featuring pair-factorized steady states (PFSS) in the closed 
system. In this model, the steady state of the particle hopping process can be tuned to yield 
properties from the zero-range process (ZRP) condensation model to those of models with spa¬ 
tially extended condensates. By varying the particle exchange rates as well as the presence of a 
global drift, we observe a phase transition from a free particle gas to a phase with condensates 
aggregated to the boundaries. While this transition is similar to previous results for the ZRP, 
we find that the mechanism is different as the presence of the boundary actually influences the 
interaction due to the non-zero interaction range. 


Condensation phenomena are observed in a broad range 
TjOf physical processes. While they are originally associ- 

i-ated with phase transitions of matter from the gas state 

to some liquid or solid state, they are also closely related 
^ to nucleation and coarsening phenomena. Examples of 
condensation appear in processes such as the formation of 
l/~) breath figures l], Bose-Einstein condensation poly- 
^mer aggregation [3], but in a wider sense also in more 
generic systems like networks as the formation of clusters 
. [2] through the accumulation of links on sites, 
o For many such systems, the involved condensation pro¬ 
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cess can be modeled as a stochastic transport process 


with a set of particles occupying a number of discrete 
• • sites. With particles representing microscopic to macro- 
. sco Pi c objects and appropriate dynamics a wide spectrum 
^ of physical processes has been studied. Examples include 
$—i refs. [2||4] mentioned above, but also processes such as 
' u wealth condensation [5] or traffic flow [6j. 

The zero-range process (ZRP) with condensation dy¬ 
namics [t}[ 8] is a well known paradigm of such transport 
processes. While it has a fully symmetric steady state, 
above some critical density p c the symmetry breaks spon¬ 


taneously and a particle condensate emerges at a single 
site such that the density at the remaining sites stays 
critical. When short-range interactions are introduced, 
a similar condensation process can be observed with the 
main difference, that condensates can be spatially ex¬ 
tended [9 12 . In this work, however, we shall consider 
condensation not as an effect in the steady state, but as 
a signature of a boundary induced phase transition. Such 
transitions can occur in driven systems, where the drive is 
implemented in terms of the interaction at the boundary 
of the system [l3]. A well known transport process with 
such a transition is the totally asymmetric simple exclu¬ 
sion process (TASEP) EES , where a high-density, a low- 


density and a maximal-current phase exist 16 . For the 


ZRP, such effects of open boundaries and driven particle 
exchange have been studied using analytical and numeri¬ 
cal methods by Levine et al. 


17 


While the interaction of particles with the boundary of 
the ZRP is merely the injection and removal of particles, 
the presence of the boundary influences hopping at nearby 
sites once short-range interactions are introduced. Within 
this paper, we use the term open boundaries in the sense 
that the boundary sites only have one interaction bond 
towards the bulk of the system. While this is a natural 
choice for an isolated system and consistent with ref. 17 
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Fig. 1: Schematic representation of the dynamics of a parti¬ 
cle hopping process on a one-dimensional lattice with L sites, 
hopping rate m and drift parameters p, q. At the boundary 
sites i = L and 1 they are replaced by the removal parameters 
Pout and q ou t, respectively. Particle injection (rate parameters 
Pin, Qin) is independent of the hopping rate. 


it is also conceivable that the boundary sites interact with 
a mean-field occupation outside the system, for example 
to model a compartment of a large system separated by 
membranes. 

We consider a one-dimensional lattice with L sites and 
a gas of indistinguishable particles. Each site i can hold 
an arbitrary number rrii of particles. Particles can hop be¬ 
tween sites as well as enter and exit the system through the 
boundaries. The dynamics in terms of a discrete stochastic 
process consists of two steps as indicated in fig. [T] First, 
a particle may leave from a randomly chosen site i with 
probability proportional to the hopping rate rq. Second, 
the particle hops to the left with rate q or to the right 
neighbor with rate p. This allows the implementation of 
symmetric (p = q = 1/2, p in = g in ,Pout = 9out) as well as 
partially and totally asymmetric [p = 1 , q = q m = q 0 u t = 
0) hopping dynamics. At the boundary sites i = L and 
1 , the drift parameters p and q are replaced by the re¬ 
moval parameters p out and g out which enter into the rates 
of particle removal ulPova and u\q ou t, respectively. The 
injection rates are independent of the occupation num¬ 
bers and hence directly given by the parameters pi n and 
q m . Due to this explicit particle injection and removal 
through the open boundaries, the total number of parti¬ 
cles M(t) = J2f=i m i is n °t conserved. For the hopping 
rate function, we consider the form 


Ui = 


n 


g(rrii - l,m,j) 
g{m i ,m j ) 


(1) 


with a symmetric, non-negative weight function g(m, n ) 
given for each bond (i,j) of the lattice^ Whereas in a 
closed system with the total number of particles conserved 
this choice leads to a steady state |9 ,[Tl of the form 

PM,L({m}) = — - TT g(mi,mj)Sj2^ , m . M > ( 2 ) 

4M,L , , ‘- 1 

<*.J> 

this is not necessarily the case with open boundaries. How¬ 
ever, since the steady state of the closed system is rooted 

1 This form of the hopping rate with fixed i implicitly realizes the 
fact, that there is no interaction beyond the boundary sites. 


Fig. 2: Regimes of different condensate shapes of the considered 
tunable model: No condensation for 7 > 1, single-site conden¬ 
sates for P < 7 < 1 , rectangular condensates for 7 < fi < 1 
and smooth parabolic condensates for ft > 1 and 7 < 1 . 


in the factorization property over bonds, i.e., pairs of sites, 
one usually still refers to it as pair-factorized steady state 
(PFSS). The normalization constant Zm,l in the steady 
state © has the same function as the partition function 
in an equilibrium model. 

The weight function g(m, n) gives the interaction be¬ 
tween particles. Here, we consider the tunable weights 


g{m, n) = exp 


| m — n\ — -( m 7 + n?) 


(3) 


proposed by Waclaw et al. 10 11 . By tuning the param¬ 


eters /3, 7 , properties of the condensation process such as 
the critical density and the condensate’s shape and ex¬ 
tension can be chosen 12 . This allows us to study the 


model in a regime with strong nearest-neighbor interac¬ 
tions, where spatially extended condensates occur in the 
steady state model (7 < 1 and 7 < /3 < 1 for rectangu¬ 
lar and f3 > 1 for smooth condensates), while still being 


17 


able to directly compare to the results of Levine et al. 
by tuning the model towards the weak nearest-neighbour 
interaction regime (/3 < 7 < 1 ), where its properties are 
similar to the ZRP considered there. Figure [2] summarizes 
the phases of this model with periodic boundaries. 

From the research in ref. 17 we know that the ZRP 


condensation model with hopping rates u(m) = 1 + b/m 
in the condensation regime b > 2 has a phase transition 
induced by boundary drive. A homogeneous gas phase 
where a steady state still exists and one (symmetric hop¬ 
ping) or two (totally asymmetric hopping) phases with 
single-site condensates aggregated to the boundary are ob¬ 
served. This transition occurs for p m > p out for symmetric 
and pi n > 1 for totally asymmetric hopping. For the ZRP, 
the presence of the boundary does not affect the inter¬ 
action at the boundary sites other than by injection or 
removal of particles. This makes it possible to understand 
the condensate formation at a boundary site for totally 
asymmetric hopping by only considering the difference of 
the fluxes into and out of that site 17 . The occupation 


number of the boundary site is then treated as a biased 
bounded random walk with drift equal to this net flux. 
As shown in fig. [3j this drift becomes positive for suffi¬ 
ciently high influx rate p- ln > 1 and occupation m\ so that 
a condensate can emerge and grow after a sufficiently large 
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Fig. 3: Drift m 2 ) of the occupation number of the 

first site, mi, for different occupations at the second site, m 2 , 
for our model (/? = 1 . 2,7 = 0 . 6 ) at p m = 1 as well as for weak 
short-range interactions (/? = 0.4 ,7 = 0.6, m 2 = 0) and the 
ZRP with totally asymmetric (p = l,q = 0) dynamics. For 
/3 = 0.4, there is no significant dependence on the second site 
occupation m 2 , for the ZRP there is none. 


fluctuation of mi. For the ZRP, the effect of these bound¬ 
ary condensates to the bulk system is then, that they act 
as particle reservoirs that hold the particle density in the 
bulk stationary at the critical density p c = l/(b— 2). It 
is also apparent from fig. [3] that with short-range inter¬ 
actions this drift behaves entirely differently and also de¬ 
pends on the second site’s occupation m 2 . We will pick up 
this observation further below, when we discuss the corre¬ 
sponding phase of aggregate condensate formation for the 
short-range interaction model. 

In our own preliminary work 


18 we additionally ob¬ 


served for the PFSS model the formation of condensates 
in the bulk system in the aggregate condensate phase for 
symmetric hopping. From this we could also assume that 
with short-range interactions the phase diagram induced 
by particle exchange and external drive is similar to that 
of the ZRP. However, there were limitations of the simu¬ 
lation method of the stochastic process because no upper 
bounds exist for the hopping rate function generated by 
the weights ( 0 ) , so that we were not yet able to study the 
phase diagram with short-range interactions to our satis¬ 
faction. 

To avoid these problems, we used here an improved 
rejection-free kinetic Monte Carlo (KMC) method, orig¬ 
inally introduced for the simulation of coupled rate equa¬ 
tions in chemical systems by Gillespie 19|20 . This allows 
us to directly use the interaction rates of the system which 
makes simulation much more efficient in situations where 
no steady state exists. 

In the following we will outline how we identify the 
phase diagram. As a main observable to indicate the ex¬ 
pected phase transition we measure the total number of 
particles M ( t ) as shown for symmetric hopping in fig. [4] 
and estimate the scaling exponent a under the assumption 
that M(t) oc t a grows as some power of time due to ab¬ 
sorption of a fraction of particles entering the system. We 
can thus identify regions in (pi n ,Pout), where the system 
continually absorbs particles and the steady state breaks 



Fig. 4: The total number of particles M(t) in the system ver¬ 
sus time. For sufficiently large times, a stationary base line is 
a signature of the gas phase while linear growth occurs in the 
aggregate condensate phase. Curves averaged from 25 inde¬ 
pendent simulations at /3 = 1.2, 7 = 0.6, p = q = 1/2, L — 256. 

down in the regions where a « 1. This is demonstrated in 
fig-! which shows the average total number of particles 
versus time for several combinations of particle exchange 
rates p in and p ou t- 

In contrast, for sufficiently small values of the influx 
rate p; n we obtain a ss 0 which reflects that the number 
of particles in the system remains stationary and suggests 
that a steady state exists. In the latter case, the particles 
are homogeneously distributed in the system as a thin gas, 
rarely interacting (phase G). This is visible in the low av¬ 
erage occupation number density pbuik hr fig. [0| as well as 
in the average occupation profile shown in fig. [T| (flat pro¬ 
files) . In the gas phase the overall density p is identical to 
the bulk density where any particle accumulations at the 
boundaries are neglected. In the regions with positive a 
(dark regions in fig. [ 5 ]), particles accumulate in the system. 
This occurs through the emergence of spatially extended 
condensates at the boundaries, as can be seen in fig. [7] As 
in this regime the condensates are bound to the bound¬ 
ary, we refer to it as the aggregate condensate phase (A). 
While for symmetric hopping (p = q = 1/2) identical for¬ 
mations appear at both boundaries [fig. [7j a)], distinctly 
shaped and independent condensates are observed for to¬ 
tally asymmetric hopping (p = l,q = 0) [fig. [7)(b)]. These 
condensates in the latter case appear separately in the re¬ 
gions Ai n and A out or together in the region A of figs. [5] 
and! 

At the boundary with particle influx, for a sufficiently 
high injection rate p; n in the regime A; n , an aggregate 
condensate forms that obeys the same envelope shape as 
a bulk condensate in the periodic model [ll [both bound¬ 
aries in fig. 0a), left boundary in (b)]. There also appears 
a regime A out , where particles condense prior to leaving 
the system, as well as a coexistence region A [for example 
Pin = l,p 0 ut = 0.2 in fig. 0b)]. For symmetric hopping 
the inbound condensate coexists with the outbound con- 
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Fig. 5: Scaling exponent a(pi n ,p out ) of the total mass M(f) 
over time for (a) symmetric (p = q = 1/2,pi n = ?i n ,Pout = 
(/out) and (b) totally asymmetric (p = 1 ,q = q m = (/out = 0 ) 
dynamics estimated at large times 10 7 < t < 10 s for /3 = 
1 . 2,7 = 0 . 6 . 


Fig. 6: Bulk density p bu ik(Pm,Pout) for (a) symmetric (p = 
q = 1 / 2 , Pin = (/i n ,Pout = (/out) and (b) totally asymmetric p = 
1 , q = gin = g ou t = 0) dynamics with interaction parameters 
P = 1 . 2,7 = 0.6. The critical density observed in a closed 
system of comparable size L = M = 256 for these parameters 
Is Pcrit SS 0.3. 


densate, the former, however, dominates the shape due to 
its much larger contribution. The transition lines from the 
gas phase G to the aggregate condensate phase A for the 
coupling constants /3 = 1 . 2,7 = 0-6 are found at p- m > 
0.625p O ut for symmetric and p m > Pm jC rit = 1-43 ± 0.02 
for totally asymmetric hopping as shown in figs. [5]and[6j 
We also observe that the location of the transition lines 
depends strongly on the couplings /?, 7 as well as the type 
of interaction at the boundaries. For instance, by decreas¬ 
ing the strength of the short-range interactions, the criti¬ 
cal lines approach those observed for the ZRP until they 
match for /3 < 7 . An interesting modification of the inter¬ 
action at the boundaries that consists of adding an explicit 
bond from the boundary sites to an outer site with occu¬ 
pation 771 oo, effectively introducing fixed boundary condi¬ 
tions. leads to a shift of the transition lines towards larger 
values of p 0 ut for asymmetric hopping, so that the phases 
A out and A become significant regions of the phase dia¬ 


gram 
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While the phase diagrams given in figs. [5] an d© are simi- 
lar to those of the ZRP 117,18 with comparable boundary 
conditions, we do observe qualitative differences. With 
short-range interactions, the bulk density of the gas phase 
increases towards the critical density while approaching 
the transition line towards aggregate condensate forma¬ 
tion (see fig. [ 6 ]). However, while the occupation density 
in the bulk system, apart from the aggregate condensates, 
remains at criticality for the ZRP as well as weak short- 
range interactions [5 < 7 , we observe a sharp decline at 
the transition line to a bulk density considerably smaller 
than the critical density here. 


On a higher level, with significant short-range interac¬ 
tions ([3]), i.e. /3 > 1, the formation of the aggregate con¬ 
densate at the boundary sites for asymmetric dynamics 
cannot simply be understood by means of a bounded, bi¬ 
ased random walk of the first site’s occupation number 
as with the ZRP. Only with weak nearest-neighbour cou¬ 
plings j3 < 7 in our considered model this explanation 
remains valid, as suggested by the monotony and limit 


of the drift shown in fig. [3] With short-range interac¬ 
tions such a drift also depends on the occupation of the 
neighbour site. As shown in fig. [3j for /? > 1 this drift 
is positive even below the observed critical particle influx 
rate p;„ < Pi n ,crit ~ 1-43 for 1 < m 1 < 7712 but becomes 
strongly negative for m 1 > m 2 - That is, no fluctuation of 
the first sites occupation alone can lead to growth of the 
aggregate condensate in this approach. However, the re¬ 
gion in which the drift remains positive expands for higher 
occupations of the second site: The aggregate condensate 
can only emerge and exist as a spatially extended conden¬ 
sate because a sufficiently high occupation at a couple of 
the first sites of the system is required to result in positive 
drift of the first site’s occupation number mi. 

A second mechanism of the emergence of the aggregate 
condensate is only observed with short-range interactions. 
For the ZRP or /? < 7 the condensate directly forms at 
the boundary site. With strong short-range interactions 
(3 > 1, however, the condensate may initially emerge in 
the bulk system and then relatively fast connect to the 
boundary. 

We observe boundary drive induced phase separation 
that becomes visible as the formation of aggregated con¬ 
densates at the system boundaries from a homogeneous 
gas phase. The phase diagram itself is similar to that of 
the ZRP as discussed in ref. 17 , specifically the proper¬ 


ties of the homogeneous gas phase where the steady state 
is not broken are much alike. The observed aggregate con¬ 
densates are spatially extended as expected. The accumu¬ 
lation of particles in these condensates is much stronger 
than in the ZRP or than expected, decreasing the particle 
density in the bulk system far below the critical density of 
the steady state model, therefore excluding the formation 
of stable droplets in the bulk system. 

Since with short-range interactions nearby boundaries 
do affect the dynamics, it is worthwhile to also look at 
the influence of the specific implementation of the bound¬ 
aries. One approach to this could be to assume virtual 
bonds across the boundaries to a mean-field occupation 
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(Pin,Pout): * (0.7, 0.2) - (0.7,1.2) • (0.7, 2.0) - (1.0, 0.2) - (1.0,1.2) • (1.0, 2.0) * (2.0, 0.2) - (2.0,1.2) • (2.0, 2.0) 


Fig. 7: Average occupation number profile for (a) symmetric hopping (p = q = 1/2,pi n = Qm,Pout = <7out) and (b) totally 
asymmetric hopping (p = l,q = qi n = g out = 0) for varying particle influx and outflux rates. The interaction parameters are 
/? = 1.2 ,7 = 0.6 and the system is of size L = 256. 


to oo as suggested above, so that the homogeneity of the 
system is not a priori broken by the missing bonds. Al¬ 
though this is a large change of the interaction strengths 
at the boundary, we merely observe different positions of 
the phase boundaries as well as a deformation of the shape 
of the aggregate condensate as it couples to the mean-field 
occupation m 00 beyond the boundary. A more interesting 
change is to make particle injection and removal symmet¬ 
ric. With the rate of particle removal set constant to p ou t, 
we observe an additional phase with a steady state for 
symmetric hopping in between the gas and aggregate con¬ 
densate phases. Most notably, it features a single large 
condensate that extends the complete bulk of the system 
but falls to zero nearby the boundaries. For a more de¬ 
tailed discussion of these observations, we refer to ref. 21 


cific interaction with the boundary strongly influences the 
observable phenomena to the point where extra phases 
become observable. 

As an outlook to future work we see many possibilities 
for further research. For example, when we considered 
constant removal at the boundary to achieve a symmetry 
between particle injection and removal, one might addi¬ 
tionally take the other choice and use hopping injection, 
i.e., particles must successfully hop into the system. An¬ 
other intriguing variation is to actually implement the sys¬ 
tems considered here embedded in a larger, possibly pe¬ 
riodic, system and interpret the boundaries of the inner 
system as membranes. 


* * * 


To conclude, we numerically studied the emergence 
of phase transitions induced by driven particle exchange 
through open boundaries in a stochastic transport process 
with tunable short-range interaction. We observed a gas 
phase and two aggregate condensate phases, where parti¬ 
cles accumulate at the boundaries. We presented phase 
diagrams for the extended smooth condensate regime of 
the model and discussed the phases’ properties. We com¬ 
pared these observations with analytical and numerical re¬ 
sults for the zero-range process by Levine et al. 17 where 
similar phases are observed. While the gas phases in both 
systems are qualitatively identical, there are noteworthy 
differences in the aggregate condensate phases such as the 
spatial extension and in the case of asymmetric hopping 
different shapes of the aggregate condensates and the be¬ 
low critical bulk density. Furthermore, we found that the 
mechanism of the formation of aggregate condensates is 
different as well as supplemented by an additional mecha¬ 
nism, where the condensate forms in the bulk and quickly 
aggregates at the boundary. Finally, we shortly discussed 
other types of open boundaries that we deem important for 
a systematic study of a short-range interaction model with 
condensation where we observed additional phases featur¬ 
ing a large bulk condensate and a homogeneous fluid, re¬ 
spectively. Therefore, additionally to the boundary drive 
that our model shares with the ZRP model 17 , the spe¬ 
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